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Abstract. We introduce the notions of d-orthostochastic, <i-unistochastic, 
and d-qustochastic matrices. These are the particular cases of ¥ d - 
bistochastic matrices for F = R, C, H. The concept is motivated by 
mathematical physics. When d = 1, we recover the orthostochas- 
tic, unistochastic, and qustochastic matrices respectively. This 
work exposes the basic properties of F d -bistochastic matrices. 
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1. Introduction 

A square matrix P with nonnegative entries adding up to one 
in every row and column is called a bistochastic matrix. Bistochastic 
matrices come up in probability, combinatorics, mathematical physics, 
geometry, optimization, etc [U [16J EJ H2] . Let F stand for M, C or H, let 
nGN, and let B n (resp. U(¥, n)) denote the set ofnxn bistochastic 
matrices (resp. the group of linear isometries of F n ). The matrix of 
squared norms of entries of V £ U (F, n) is bistochastic, yielding the 
squared norm mapping v : U(¥,n) — > B n . 

The ranges of v for F = R, C, H are the sets of ortho stochastic, 
unistochastic, qustochastic n x n matrices, denoted by Q n ,U n ,$) n re- 
spectively. The obvious inclusions D n C U n C Sj n C B n are proper 
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if n > 4 [5]. There are connections between these matrices and (gen- 
eralized) numerical ranges and numerical shadows [2U US EH El El 
[91 [TOl HH [201 [IH HI- They (especially unistochastic matrices) are of 
importance in quantum physics |H El [20], [8] . There is a considerable lit- 
erature regarding orthostochastic, unistochastic, and qustochastic ma- 
trices [HI El [21 [8] . However, several basic questions remain open [3 [8]. 

The motivation for the present work comes from mathematical physics. 
The entries of a bistochastic matrix P G B n are the probabilities of 
transition between the states of a classical physical system. Transi- 
tions amplitudes in quantum physics are complex numbers, and the 
transition probabilities are their squared absolute values. The tran- 
sitions amplitudes of a quantum system with n basic states form a 
unitary n x n matrix U. Let P = v(U) be the matrix of squared 
norms of the entries in U. Then P is a bistochastic matrix describing 
the corresponding classical system. Thus, a bistochastic matrix P is 
unistochastic if the classical system described by P can be quantized, 
yielding a quantum physical system whose transition amplitudes are 
the entries of a unitary matrix U satisfying v{U) = P. 

Suppose now that the basic states of a quantum system have inter- 
nal degrees of freedom. Assume, for simplicity, that the number, say 
d > 1, of internal degrees of freedom is the same for all states. The 
transition amplitudes become vectors v\ G C . This quantum system 
corresponds to a n x n matrix V = [vj] with vector entries. The uni- 
tarity condition says that the operator V : C n — > C nd is an isometry. 
Let Iso(C, n, d) denote the set of these isometries. For V G Iso(C, n, d) 
the corresponding bistochastic n x n matrix P is given by p\ = \ \vl\\ 2 . 
Thus, our quantum systems correspond to V G Iso(C, n,d); the reduc- 
tion from a quantum to the classical system is given by the squared 
norm map v : Iso(C, n, d) — > B n . Its range is the set OS(C, n, d) C B n 
of d-unistochastic matrices. A classical system described by P G B n 
admits a quantization with d internal degrees of freedom if and only if 
P G OS{C, n,d). 

The same construction, based on R (resp. H) leads to d- orthostochastic 
matrices OS(R,n,d) (resp. d- qustochastic matrices OS(M, n, d)). To 
streamline the exposition, we will refer to OS(¥, n, d) C B n as the set 
of (F, n, d) -bistochastic or simply W d -bistochastic matrices. 

We have structured the exposition as follows. In section [2] we pre- 
cisely define the above notions and expose a few basic properties of 
vector F d -bistochastic matrices. Section E] continues this exposition. 
Section H]is devoted to (n— l)-orthostochasic matrices. It contains our 
main results, Theorem [1] and Corollary [21 and a few examples. In the 
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concluding section [5] we state a conjecture regarding (i-orthostochastic 
matrices. 



2. Definitions and basic properties 

A m x n matrix is a collection of mn entries organized in m rows 
and n columns. In general, we will denote matrices by capital letters, 
and denote their entries by the corresponding low case letters with 
subscripts and superscripts. For instance, let P = [pi] be a n x n 
matrix. The lower (resp. upper) index stands for the column (resp. 
row). The entries of a matrix may be numbers, vectors, etc. 

We will denote by F any of the fields 1 or C or H. Recall that a real 
n x n matrix P = \pl] is bistochastic iip\ > and Y^j=i Pi = Yl^i Pi = 
1. We will denote by B n the set of bistochastic n x n matrices. 

For N G N let F^ denote the vector space of iV-tuples over F. Let 
d,neN. The space ¥ nd , decomposed as the direct sum of n copies of ¥ d 
will be denoted either by ©" =1 F^ or simply by ¥ nd if the decomposition 
is clear from the context. 

Let M(¥ d ,n) be the space ofnxn matrices with entries in F . 

Lemma 1. A matrix V = [vj] G M. (¥ d , n) satisfies for all 1 < j, k < n 
the equation 

n 

(1) z>i,«*> = *0'.*)- 

i=l 

if and only if for all 1 < j, k < n we have 

n 

(2) D^ v ?> = *o-,*)- 

i=i 

Proof. Let ( , ) n and (,) n d denote the scalar products in F n and ¥ nd 
respectively. Any V G M. (¥ d , n) determines a linear operator V : F n — > 
¥ nd . The adjoint operator V* : ¥ nd ->- F n satisfies for all x G F ra , y G ¥ nd 
the identity 

(3) (Vx, y) nd = { x, V*y) n . 

Let Idfc denote the identity operator on ¥ k . Equation ([1]) means that 
V : F n — >■ ¥ nd is an isometry, i.e. 

(Vx,Vy) nd = (x,y) n . 



By equation ()3]), this is equivalent to V*V = Id n , i.e., equation ()2]). □ 
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Let Iso(F, n, d) denote the set of matrices V G Ai(¥ d , n) satisfying 
the equivalent equations ([I]), (T5]). For V = [v?] G Iso(F, n, d) we define 
the real n x n matrix P by 

(4) rf = (4 A) = IKII 2 - 

By equations (OQ) and J|, Pe S n . Thus, equation fl4]) defines a map- 
ping 

z/ : Iso(F, n, d) — >• i3 n . 

Definition 1. .4 matrix P G £> n zs (F, n, <i)-bistochastic (¥ d -bistochastic 
for brevity) if P = u(V) for some V G Iso(¥, n, d). 

When d = 1 and F = R, C, EI Definition [T] yields orthostochastic, 
unistochastic, and qustochastic n x n matrices respectively. We will 
denote by OS(¥, n, d) C S n the set of F d -bistochastic matrices. 

Proposition 1. 1. The set OS(W,n,d) C B n is closed. 

2. There are inclusions 

OS(R,n,d) C OS(C,n,d) c C5(H, n, d), 05(F, n, d) c 05(F,n,d+l). 

3. We have 

OS(¥,n,n) = B n . 

Proof. The set Iso(F, n, d) is a compact manifold and v : Iso(F, n, d) — > 
B n is a differentiable map, yielding claim 1. Claim 2 is immediate from 
the definitions. We will prove claim 3. Let P = [p^] G B n and let d G N 
be arbitrary. Then P = u(V),V G Iso(F, n, d) if there are vectors 
v\ G ¥ d such that \\vj \\ 2 = pj and the vectors 

w i = (vl,...,v?)e¥ nd , l<i<n, 

are pairwise orthogonal. We have 

(5) (wi, Wj) nc i = (fj 1 , vj) n H h «' , v") „ . 

If for 1 < k < n the vectors vf, v%, . . . , v% are pairwise orthogonal in ¥ d , 
then, by equation (jSJ), the vectors . . . , w n are pairwise orthogonal in 
¥ nd . When d > n, the space ¥ d contains n pairwise orthogonal vectors 
with arbitrary norms. □ 

The following is immediate from Proposition [TJ 

Corollary 1. Let ¥ be any o/M, C, EL Then the following holds. 

1. For any n G N t/iere zs a unique <i m i n (F, n) G N stzc/i £/ia£ 0«S(F, n, (i m in) 
B n and OS(¥, n, d) ^ B n for d < d min . 

2. We have 

d m in (EI, n) < d min (C, n) < (i min (M, n), d min (F, n) < d min (¥, n + 1), 



GENERALIZED UNISTOCHASTIC MATRICES 5 

and d m i n (F, n) < n. 

Let now P G B n . By Proposition CEJ there exist d < n such that 
P G OS(¥,n,d). Let d m i a (P, F) be the minimal such cf. Then 

(6) d min (F, n) = max{d min (P, F) : P G 

We will informally refer to d m i a (P, F), d m in(F, n)) as the minimal num- 
ber of internal degrees of freedom for F- quantization. 

There are obvious identifications: 

Iso(R,n,l) = 0{n), Iso(C,n,l) = U(n), Iso(H,n, 1) = Sp(n). 

and 

OS(R, n, 1) = O n , OS{C, n, 1) = il„, OS{M, n, 1) = S) n . 
The following is well known: 

OS(R, 2, 1) = OS(C, 2, 1) = OS(M, 2, 1) = B 2 , 

OS(R, 3, 1) c C5(C, 3, 1) = OS(W, 3, 1) c B 3 
and the inclusions are proper. For n > 3 there are proper inclusions [5] 

05(R, n, 1) C OS(C, n, 1) c OS(U, n, 1) c B n - 

3. Relationships, symmetries, and dimension count 

There are several relationships between (F, n, <i)-orthostochastic ma- 
trices for various values of F and d. 

Proposition 2. For any n > 1 and d G N there are natural inclusions 
OS{C,n,d) C 05(R,n,2d), 05(H,n,d) c C5(C,n,2rf) 

and 

OS(m,n,d) C 05(R,n,4d), 

Proof. The decomposition z = x + \J—\y identifies C d and M 2d . Let 
( , )c and ( , )r be the complex and the real scalar product on C d re- 
spectively. The relationship 

$t((u,v) c ) = (u,v)jh 

yields the proper inclusion Iso(C, n, d) C Iso(R, n, 2d). It is compatible 
with the squared norm maps v : Iso(R, n, 2d) — > B n , v : Iso(C, n, d) — > 
B n , yielding the first inclusion. The second follows similarly from the 
isomorphism M d = C 2d , and the third inclusion is the composition of 
the former two. □ 
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Let W n C M. (R, n) be the group of permutation matrices. It acts 
on Ai (R, n) by left and by right multiplication, yielding the action 
of W n x W n on M(R,n) which preserves i3 n . Let Iso(F fc ) be the 
group of linear isometries of ¥ k . The set Iso(F, n, d) is invariant un- 
der precomposition (resp. postcomposition) with elements in Iso(F n ) 
(resp. Iso(F na! )) yielding the action of Iso(F" d ) x Iso(F n ) on Iso(F, n, d). 
Let Diag(F n ) C Iso(F n ) (resp. Diag(F,n,d) C Iso(F nd )) be the sub- 
group of diagonal (resp. block-diagonal) isometries. By restriction, 
the subgroup Diag(F, n, d) x Diag(F n ) acts on Iso(F, n, d). The group 
W n embeds naturally into Iso(F n ) (resp. Iso(F n<i )), by permutation 
(resp. block permutation) matrices yielding an action of W n x W„ on 
Iso(F, n,d). 

Proposition 3. The squared norm map v : Iso(¥, n, d) — > OS(¥,n,d) 
is invariant under the action of Diag(¥, n, d) x Diag(¥ n ) and equivariant 
for the actions ofW n x W„ on Iso(¥,n,d) and OS(¥,n,d). 

Proof. Let V = [vj] G Iso(F, n, d). The action of Diag(F n ) multiplies 
the n vectors v\ 6 ¥ d , where i is fixed and 1 < j < n, on the right by 
the same element A, G F with ||Aj|| = 1. The action of Diag(F, n, d) 
multiplies the n vectors v\ G ¥ d , with j fixed and 1 < i < n, on the 
left by the same isometry Uj G Iso(F d ). These actions do not change 
the norms of vectors v\. 

If P = v(V) thenp^ = 1| 2 . The action of W n x W n on Iso(F,n, d) 
permutes the rows and the columns of V the same way as its action on 
B n permutes the rows and the columns of P. □ 

By the preceding discussion, the group Diag(F, n, d) x Diag(F n ) nat- 
urally acts on Iso(F, n, d). We denote by VC(¥, n, d) the quotient space, 
i.e., VC{¥,n,d) = Diag(F, n, d)\Iso(F, n, rf)/Diag(F n ). By Proposi- 
tion 131 the squared norm map v : Iso(F, n, d) — > B n uniquely descends 
to a mapping of T>C(F, n, d) which we will also denote by v. Thus, 
OS (F, n,d) = v (VC(¥, n,d)) C B n . Let dimX denote the real dimen- 



sion. 



Proposition 4. The following equations hold: 




d 2 - 



d + 1 



2 



n. 
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dim (VC(C, n, d)) = (2d - l)n 2 - (d 2 + l)n + 1, 



and 



(9) dim (VC{W, n, d)) = (Ad - 2)n 2 - (d 2 + d + 2)n. 
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Proof. We have 



dim (Iso(R, n, d)) = (d - -)n 2 - -n, dim (Iso(C, n, d)) = (2d - l)n 2 , 



dim (Iso(H, n, d)) = {Ad - 2)n 2 + n. 



Specializing to d = 1, we recover the well known formulas 



dimO(A;) = — J -,dimU(k) = k 2 , dim Sp(k) = 2k 2 + k. 
The groups Diag(F") satisfy 




Diag(M n ) = {±l} n ,Diag(C") = U(l) n , Diag(H") = Sp(l) n , 



The actions of Diag(F, n, d) and Diag(F n ) on Iso(F, n, d) are free and 
commute. They are transversal, except that both Diag(C, n, d) and 
Diag(C n ) contain U(l) as the group of scalar unitary matrices. This 
information and the above formulas yield the claims. We leave details 



4. On (n - 1)-orthostochasticity of n x n matrices 

Let P e B n . By Proposition (TJ d min (P,R) < n. We will show 
that if P satisfies mild non-degeneracy assumptions and n is odd, then 
d m in(P, K.) < n — 1. Let n e N. In what follows we will use the cyclic 
convention for indices: i + n = i. 

Lemma 2. Let n G N be an odd integer. Let £i, . . . ,£ n be arbitrary. 
Then the system 

(10) Xi + x i+1 = Ci+2, 1 < i < n, 
has a unique solution 

(11) 2x { = + + ^+2 - ^+3 + ^+4 — I h €i+n-i, 1 < i < n. 

Proof. The matrix of the linear system ( ITOj) is nondegenerate, hence 
the solution is unique. The reader will easily verify that equation (iTTj) 
yields a solution. □ 

Lemma 3. Let n e N be an odd integer. Let ^i,...,^ n be positive 
numbers satisfying the inequalities 

ii + 6+3 + 6+5 + • • ■ + £i+n-2 < 6+1 + 6+2 + 6+4 + " " " + 6+n-l 

for 1 < i < n. Then there exists a real skew- symmetric matrix A = [aj] 
satisfying 



to the reader. 



□ 



(12) 




Ki<n 
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Proof. Set 0? = (a]) 2 . Then B = is a symmetric n x n matrix with 
non-negative entries. The matrix A satisfies equation ( TT2l) if and only 
if B satisfies 

(13) Yl % = £i, 1 < 3 < n. 

\<i<n 

For xi, . . . , x n G M. + set 

b\ = b\ = x n , b\ = b\ = xi, 6| = 6^ = x 2 , . . . , b n ~ l = b n n _ x = x n _ 2 , 

Hi = ^i = x n-i, and let 6^ = for all other pairs of indices. Then B 
satisfies equation (TT3"j) if and only if x%, . . . ,x n satisfy equation ( ITU]) . 
The claim now follows from Lemma HJ □ 

Theorem 1. Let n > 1 be an odd integer. Let P = [p^] be a n x n 
ortho stochastic matrix. Suppose that P satisfies the inequalities Pij^O 
for i j and that for 1 < i < n we hav^ 

(14) Pt+Piti+Pts + • • • +p£K=! < pttl+pitl+Pn-t + ■ ■ • +*ffi=i- 

Then P is (n — 1)- ortho stochastic. 

Proof. We will find V = [vj] G Iso(M, n, n — 1) such that P = v{V). For 
1 < j < n let v[, . . . , Vj_ X i v j+n ■ ■ ■ -> v h e K n_1 be mutually orthogonal 
vectors such that \\p 3 k \\ 2 = p{- The vectors v\ = ^(p*.) -1 ' 2 for 1 < k < 
n,k 7^ j, form an orthonormal basis in M n_1 . 

For 1 < i < n let Wi G M"^™ -1 ) and the notation ( , )„(„_!) be as in 
the proof of Proposition [TJ For 1 < j < n we set 



(15) vi = J2 



Thus, Wi, . . . ,w n are the column vectors of V = [vf]. The equation 
v{V) = P is equivalent to the two systems of quadratic equations: 

(16) (w h lOj)n(n-l) = 
for i 7^ j and 

(17) \\v]\?=P i j 

for 1 < j < n. 

By equation (TT5l) . we have 

(18) «;j) n („_i) = a • + a) 



1 We use the cyclic convention for indices. 
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for % ^ j and 



(19) ikh 2 = e^) 2 - 

t=i 

Thus, it suffices to find a skew-symmetric n x n matrix A = [a?] 
satisfying 

n 

(20) »:s 2 4 

Lemma [3] yields a solution. □ 

Remark 1. The proof of Theorem [1] yields an explicit solution of the 
equation u{y) = P. Let ei,...,e n _i G M n_1 be the standard or- 
thonormal basis. For 1 < j < n set v 3 k = (j^e) e k if k < j and 
v l = (Pk) l ^ e k-i ^ k > 3- The numbers a\ defining the vectors 
v J j G M 71-1 satisfy = if i — j ^ ±1 mod n. For pairs i, j satis- 
fying i — j = ±1 mod n, Lemmas [2, [3] explicitly yield (a]) 2 as linear 



combinations of p k where 1 < k < n. Set a\ = \J (aj) 2 if j < i and 



(V; 



3 



il) 2 if i < j. 



Corollary 2. Let n > 1 fre an odd integer. Then 

1. T/ie sei OS{R, n, n — 1) contains all bistochastic matrices P = [pi] 
such the numbers p\, . . . ,p™ satisfy the inequalities in equation (TH1) : 

2. For an open set of P G B n we have d min (P, M) < n— 1. 

Proof. Let P = [p^] beanxn bistochastic matrix. By Theorem [IJ P G 
C«S(M, n, n — 1) if the numbers p\, . . . ,p™ satisfy the strict inequalities 
in equation ffTil) and p\ ^ for i ^ j. This is a nonempty open set, 
hence claim 2. By Proposition [TJ OS(W,n,n — 1) contains its closure, 
yielding claim 1. □ 

Example 1. Let n = 3. The open set in Corollary [2] is the set of 
P G £>3 such that 7^ for i 7^ j and 

(21) p\ <pI+pI,P 2 2 <P 3 3 +pI,pI <p\+pl 

Every 3x3 matrix satisfying these inequalities is 2-orthostochastic. 
Note that the set of orthostochastic 3x3 matrices has positive codi- 
mension in $ 3 [2J 0, H5] . Inequalities (I2T1) hold if and only if p\,p%p\ 
are the side lengths of a nondegenerate triangle. We point out that the 
triangle inequalities come up as conditions of unistochasticity for 3 x 3 
matrices [8]. 



10 



EUGENE GUTKIN 



Example 2. Let n = 5. The set of 4-orthostochastic 5x5 matrices in 
Corollary [2] consists of bistochastic matrices satisfying the following: 

pI+pI <pI + pI+pIpI + pI <pl+pt+p{, 

and 

pI + p\ < p\ + pI + pI, pI + p\ < pI + p\ + pI pI+p\<p\+pI + pI 

(These inequalities do not have an immediate geometric interpreta- 
tion.) The interior of the set of 4-orthostochastic 5x5 matrices con- 
tains [pi] € £> 5 such that pi ^ for i ^ j and the above inequalities 
are strict. 

5. Concluding discussion 

To describe precisely the sets of orthostochastic, unistochastic, and 
qustochastic n x n matrices for arbitrary n seems a difficult problem 
ref. The concept of F d -bistochastic introduced here replaces this prob- 
lem with another, seemingly simpler, but still a meaningful problem: 
To characterize F d -bistochastic n x n matrices for arbitrary n. Our 
results suggest that this question is still nontrivial but less subtle. 
The two questions are related. Note that, by Proposition [2j the set 
of 2-orthostochastic matrices contains the set of unistochastic matri- 
ces, and the set of 2-unistochastic matrices contains the set of qus- 
tochastic matrices. Let n be arbitrary, and let us vary d G N. By 
Proposition the set of F d -bistochastic coincides with B n as d reaches 
rf m in(F, n) < n. The above results and dimensional considerations sug- 
gest that d min (W,n) does not indefinitely increase, as n goes to oo. In 
particular, for F = R, the following conjecture is plausible. 

Conjecture 1. Denote by e£ m j n (n) the smallest value of d such that the 
set of d- orthostochastic matrices coincides with B n . Then there exists 
n such that for n > n we have 

2 < d mhl (n) < 3. 

The material in section 0] suggests that value of <i min (n) may depend 
on the parity of n. The author believes that the threshold dimension 
should not be too large. Most likely, no = 3. 
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